The information about the state of a charge qubit gained by a weakly coupled 

quantum point contact 
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We analyze the information that one can learn about the state of a quantum two-level system, 
i.e. a qubit, when probed weakly by a nearby detector. We consider the general case where the qubit 
Hamiltonian and the qubit's operator probed by the detector do not commute. Because the qubit's 
state keeps evolving while being probed and the measurement data is mixed with a detector-related 
background noise, one might expect the detector to fail in this case. We show, however, that under 
suitable conditions and by proper analysis of the measurement data useful information about the 
initial state of the qubit can be extracted. Our approach complements the usual master-equation and 
quantum-trajectory approaches, which describe the evolution of the qubit's quantum state during 
the measurement process but do not keep track of the acquired measurement information. 



I. INTRODUCTION 



Solid-state systems are among the most promising can- 
didates for future quantum information processing de- 
vices. One type of such systems are superconductor- 
and semiconductor-based charge qubits ^ . These qubits 
are commonly measured by capacitively coupling them to 
quantum point contacts (QPC) or single-electron transis- 
tors (SET), such that the current in the detector is sen- 
sitive to the charge state of the qubit 2, 3, 4, 5, 
By measuring the current passing through the detector, 
one can infer the state of the qubit. One limitation that 
arises in practical situations is that, in order to minimize 
the effects of the detector on the qubit before the mea- 
surement, the qubit-detector coupling strength is set to a 
value that is small compared to the qubit's energy scale 
[|. As a result one must deal with some form of weak- 
measurement regime. This type of weak, charge-sensitive 
readout works well when the qubit is biased such that 
the charge states are eigenstates of the Hamiltonian and 
therefore do not mix during the measurement. In this 
case one can allow the detector to probe the qubit for as 
long as is needed to obtain a high signal-to-noise ratio, 
without having to worry about any coherent qubit dy- 
namics (Note that, since we are mainly interested in the 
measurement process, we ignore any additional qubit de- 
coherence mechanisms in the system, which could impose 
constraints on the allowed measurement time). 

In contrast to the simple situation described above, 
when the detector weakly probes the charge state of the 
qubit while the Hamiltonian induces mixing dynamics be- 
tween the different charge states, it becomes unclear how 
to interpret a given measurement signal. Since the sig- 
nal typically contains a large amount of detector-related 
noise and the measurement unavoidably destroys the co- 
herence present in the qubit's initial state, it might seem 
that this type of measurement cannot be used to deter- 
mine the initial state of the qubit, i.e. at the time that 
the experimenter decides to perform the measurement. 
Indeed, there have been a number of studies analyzing 



the measurement-induced decoherence and the evolution 
of the qubit's state in this situation 0, H, 0, d, Q , but not 
the question of how to take the measurement data and 
extract from it information about the initial state of the 
qubit. This question is a key issue for qubit-state readout 
and is the main subject of this paper. We shall show be- 
low that high-fidelity information can be extracted from 
the measurement data, provided that additional decoher- 
ence mechanisms are weak and the readout signal can be 
monitored at a sufficiently short timescale. It turns out 
that not only the measurement result, but also the mea- 
surement basis is determined stochastically in this case. 
In spite of the uncontroUability of the measurement ba- 
sis, the measurement results correspond properly to the 
initial state of the qubit. As an example, we show how 
they can be used to perform quantum state tomography 
on the qubit. These results show that under suitable con- 
ditions and by proper analysis of the measurement data 
useful information about the state of the qubit can be 
obtained. 



II. MODEL 

We consider a system composed of a charge qubit ca- 
pacitively coupled to a QPC, as illustrated in Fig. 1. The 
qubit can be viewed as a system where a charged parti- 
cle is trapped in a double well potential and can occupy, 
and tunnel between, the localized ground states of the 
two wells. We shall denote these states by \L) and \R). 

During the measurement a voltage is applied to the 
QPC, and a current flows through it. We assume that the 
QPC measures the qubit in the {\L) , |i?)} basis; the cur- 
rent through the QPC depends on whether the qubit is 
in the state \L) or \B). We further assume that the QPC 
does not induce any decoherence in the qubit's state ex- 
cept that associated with the measurement-induced pro- 
jection. For purposes of fully characterizing the operation 
of the QPC as a detector for the state of the qubit, it is 
convenient to make a few statements about the QPC's 
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FIG. 1: (color online) Schematic diagrams of (a) a charge 
qubit measured by a quantum point contact (QPC) and (b) 
the probability distributions of the possible QPC current val- 
ues for the two charge states of the qubit. The finite widths of 
the probability distributions are a result of the finite measure- 
ment time. When the distance between the two center points 
in Fig. 1(b) \Ir — Il \ is much smaller than the widths of the 
distributions, the QPC performs a weak measurement on the 
qubit in the short interval under consideration. In plotting 
(b) we have assumed that Ir > II, which would be the case 
if the qubit is defined by an extra positive charge (e.g. a hole) 
tunneling between the two wells. Taking the opposite case, 
i.e. Il > Ir, would not have any effect on the analysis of this 
paper. 

operation when the qubit Hamihonian is diagonal in the 
charge basis and the qubit is initialized in one of the 
states of the charge basis. In this case, there is no mech- 
anism by which the states \L) and \R) mix during the 
system dynamics. If the qubit is initially in the state 
the long-time-averaged QPC current is given hy II, 
and the qubit remains in the state \L). A similar state- 
ment applies to the state \R) of the qubit, with corre- 
sponding QPC current In. The QPC current therefore 
serves as an indicator of the qubit 's state in the charge 
basis {\L) ,\R)}, as long as the qubit Hamiltonian does 
not mix the states of this basis. 

On any finite timescale, there will be fluctuations in 
the QPC current, and the observed value might devi- 
ate from II ot Ifi. The longer the period over which 
the averaging is made, the smaller the fluctuations. One 
can therefore define a measurement timescale that deter- 
mines how long one needs to wait in order to distinguish 
between the states \L) and \R). The relation between 
this timescale and the qubit 's Hamiltonian- induced pre- 
cession period separates two measurement regimes: fast 
versus slow measurement, or alternatively strong versus 
weak qubit-detector coupling. Note that this distinction 
is irrelevant when the qubit Hamiltonian is diagonal in 
the charge basis, since there is no mixing between the 
states \L) and \R) in this case. 

For the remainder of this paper, we analyze the gen- 
eral case where the qubit Hamiltonian is not necessarily 
diagonal in the charge basis. We shall use the basis in 
which the qubit Hamiltonian is diagonal, thus 

iJq - -Ea,/2, (1) 

where E is the energy splitting between the qubit 's two 
energy levels, and is the z-axis Pauli matrix. We shall 
denote the ground and excited states of the Hamiltonian 



by |0) and respectively. The states of the charge basis 
can be expressed as 

|i?)=cos^|0)+sin^|l) 

|L)=sin^|0)-cos^|l), (2) 

where f3 represents the angle between the charge basis 
and the energy eigenbasis. 

III. MEASUREMENT- AND 
HAMILTONIAN-INDUCED DYNAMICS 

We start our analysis by considering a short time inter- 
val between times t and t+St. We assume that during this 
time interval a large number of electrons tunnel through 
the QPC, such that it is natural to define a QPC cur- 
rent I{t) during this short interval. We also assume that 
a weak-measurement regime exists for a properly chosen 
value of 5t, which means that the QPC-current probabil- 
ity densities (for the states \L) and \R)) are broad and 
almost completely overlap, as shown in Fig. 1(b). For 
definiteness we shall take these probability densities to 
be time- independent, Gaussian functions. The interval 
deltat is taken to be much longer than the coherence 
time of the QPC, such that the QPC's operation during 
this interval is independent of the QPC's state at earlier 
times. Finally, we take deltat to be much shorter than 
the precession period of the qubit. 

We now construct matrices (or propagators) that de- 
scribe the qubit-state evolution depending on the ob- 
served QPC current I{t): when a given value of I{t) is 
observed in the QPC, the quantum state of the qubit is 
projected (possibly partially) according to the observed 
value. Neglecting decoherence that is not associated with 
the measurement, the projection of the qubit's state is de- 
scribed by a 2 X 2 matrix that we shall call Um [I, SI, St] , 
where SI defines the size of a finite interval of QPC cur- 
rents that we identify with a single value: 

Pq(i + St) cx Um[I, si, St] p^{t) Ul[l, SI, St], (3) 

where pq denotes the qubit's density matrix, and f rep- 
resents the complex conjugate transpose. One could say 
that with the introduction of SI we are turning the prob- 
ability distributions in Fig. 1(b) into histograms with dis- 
crete possible outcomes (This discretization will also be 
used in our numerical calculations below). In order to 
identify the appropriate form for Um [I, SI, St] , we note 
that the probability of obtaining the corresponding out- 
come is 

P[I, SI, St] - Tr{;7l,[/, SI, St]UM[I, SI, St]pq{t)}. (4) 

Let us denote by Pj[I,SI,St] the probability that the 
value / (up to the dicretization parameters SI and St) 
of the QPC current is observed given that the qubit is 
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in state j. We now find that the simplest, and in some 
sense ideal, matrix UM[I,SI,5t] that satisfies Eq. ^ is 
given by 



UM[I,6I,6t] = ^PL[I,6I,St]\L) {L\- 
^PR[I,SI,St]\R) 



(5) 



This matrix could be followed by a unitary transforma- 
tion that does not affect Eq. (j4|). Any such transfor- 
mation can be incorporated into the analysis straightfor- 
wardly. 

In addition to the measurement-induced evolution de- 
scribed by Eq. ([5|) , the qubit Hamiltonian induces a uni- 
tary evolution in the qubit 's state during the time interval 
t to t + St: taking h = 1 



UniSt] 



exp 



1 + i — Sta, 
2 



(6) 



The matrices Um[I, SI, St] and UH[St] can now be com- 
bined to give the total evolution matrix 

U[I{t),SI,St] = UM[Iit),SI,St]xUH[St]. (7) 

Note that both UM[I,SI,St] and UniSt] are approx- 
imately proportional to the unit matrix in the limit 
St 0, with lowest-order corrections of order St. The 
operators UM[I,SI,St] and UniSt] therefore commute to 
first order in St. 

When a given QPC output signal I{t) [from the initial 
time t = until t = tf] is observed, one can take the 
corresponding short-time evolution matrices explained 
above and use them to construct the total evolution ma- 
trix C/Totai[-^(i : — > t/), SI, St]. Using the unit matrix as 
the total evolution matrix for t = 0, we find that 



UTot.i[I{t:0^tf),SI,St] 



U[I{tf ^ St),SI,St] X ••• X 
U[I{0),SI,St]. (8) 



Once the 2x2 matrix UTotai[I{t : ^ tf),SI,St] is cal- 
culated, one can divide it into two parts, a measurement 
matrix UMeas[Iit : — > tf),SI,St] followed by a unitary 
transformation U-[iot[I{t ■ — > tf),5I,St]: 

UToUl{t:0^tf),SI,St] = 
URotiHt : ^ tf), SI, St] X UucsiHt : ^ tj), SI, ^ 



The matrix Umcss has the form 



(10) 



where and ]ip2) are two orthogonal states. The states 
IV'i) and 1-02) represent the measurement basis that cor- 
responds to the output signal I(t), and the parameters 
Pi are the probabilities that the outcome defined by I{t), 
SI and St is obtained given that the qubit was initially 
in the state With a simple calculation, one can see 
that the measurement fidelity is given by (see e.g. p^ ) 



Fidelity = 



Pi -P2 



Pi -I-P2 



(11) 



To summarize, the QPC output-current signal can be 
used to derive the matrix UTota.i[I{t ■ tf),SI,St]. 
This matrix can then be used to determine the measure- 
ment basis, the measurement result (i.e. ±1 along the 
measurement axis), the fidelity (or in other words the 
degree of certainty about the obtained measurement re- 
sult) and the final state of the qubit (given by the mea- 
surement result transformed by the unitary, i.e. rotation, 
matrix U-ajot[I{t ■ tf),SI,St]). Note that when the 
measurement fidelity approaches one, the final state is a 
pure state that can be determined even without any a 
priori knowledge about the initial state. 



IV. RESULTS AND DISCUSSION 

We now present the results of our numerical calcula- 
tions. The calculations were performed by analyzing a 
sequence of discrete events, with each event representing 
a time steps of size St. We also discretize the possible 
values of QPC current. We have verified that the results 
presented below are insensitive to the exact discretiza- 
tion parameters, as long as we take ESt <C 1 and there 
are a large number of possible QPC current values. The 
qubit is initialized in a given state that depends on the 
specific calculation. In each time step, the value of the 
QPC current is determined stochastically using Eq. 
Based on the obtained value, the qubit's state evolves as 
described in Eq. ([3]). Following this weak- measurement 
step, a unitary transformation of the form of Eq. ^ is 
applied to the qubit's state. After a sufficiently long QPC 
signal is obtained, this signal (in all its details) is taken 
and used to extract the measurement matrix C/ivioas ex- 
plained above. This matrix is then used to extract the 
measurement basis and fidelity. 

The strength of the qubit-QPC coupling is determined 
by the relation between two parameters in the numerical 
calculations: (1) the width, or standard deviation a, of 
the QPC-current distribution functions and (2) the dis- 
tance between the average values of these distribution 
functions {Ir — II). It is more convenient, however, to 
present the results in terms of a different parameter that 
characterizes the qubit-QPC coupling strength, namely 
i?T„j/(27r), where is the timescale needed to obtain 
sufficient QPC signal to read out the state of the qubit 
(for the time being one can think of this definition as 
applying to the case when /3 = 0; but see below). If 
one considers N of the small steps considered above, the 
standard deviation of the QPC's averaged signal scales 
as a/y/N (note that a is the standard deviation for one 
step). The measurement time Tm can now be naturally 
defined as the product of the time step St and the value 
of N at which 2a/y/N = ]Ir — Il\- The measurement 
time is therefore given by 



Aa^St 



(12) 



First, in Fig. 2 we show the spherical coordinates of 
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FIG. 2: (color online) The spherical coordinates 6 and <^ 
that define the stochastically determined measurement bases 
obtained in simulations of the experiment under considera- 
tion [note that {9, 4>) = (0, 0) for the energy eigenbasis and 
{d,4>) = (/9, 0) for the charge basis]. Each figure contains 200 
points. In Fig. 2(a) _Erm/(27r) — 0.01, i.e. deep in the strong- 
coupling regime. In Fig. 2(b) Erm/ {2tt) = 0.2 [intermediate- 
coupling regime], and in Fig. 2(c) ETm/i'^ir) = 5 [weak- 
coupling regime]. In all the figures, f3 = 7r/4. In each set, 
the figure on the left is generated using the initial state \L), 
and the one on the right is generated using the initial state 
|0). Each set is an identical pair, up to statistical fluctua- 
tions, demonstrating that the initial state plays no role in 
determining the measurement basis. 



the (stochastically determined) measurement bases for 
different levels of qubit-detector coupling. In the strong- 
coupling, fast-measurement regime (Fig. 2a), the mea- 
surement basis is always the charge basis, which is the 
natural measurement basis for the detector under con- 
sideration. As the qubit-detector coupling strength is 
reduced (Fig. 2b), the measurement bases start to de- 
viate from the charge basis, and they develop some 
statistical spread. This region could be called the 
intermediate-coupling regime. In the weak-coupling, 
slow- measurement regime (Fig. 2c), the measurement 
bases are spread over all the possible directions. 

In Fig. 3, we plot the measurement fidelity as a func- 
tion of measurement duration for three different values 
of f3, keeping all other parameters fixed. We can see that 
the fidelity approaches one for long enough measurement 
duration, regardless of the angle f3. In fact, the fidelity 
is almost independent of (3. This result shows that even 
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FIG. 3: (color online) The measurement fidelity as a func- 
tion of measurement duration for three different values of the 
angle /3 between the charge basis and the energy eigenbasis: 
P — (red; lowest line), 7r/4 (green) and 7r/2 (blue). Here 
ETm/{2n) — 5, i.e. deep in the weak-coupling regime. The 
fldelity increases as the measurement duration increases, but 
the fidelity is essentially independent of the angle /3. 



though more complicated analysis is needed to extract 
useful measurement information when /3 ^ 0, the infor- 
mation acquisition rate is only slightly affected by the 
coherent, Hamiltonian-induced precession. 

The fact that the measurement basis is generally un- 
predictable, and therefore uncontrollable, is a rather 
strange phenomenon from a fundamental point of view. 
From a practical point of view, one can wonder whether 
anything useful can be done with such measurements that 
are performed in a stochastically determined basis. If one 
absolutely requires a measurement in a given basis, mea- 
surement results in different bases would be less useful. 
One could then treat the deviation of the observed mea- 
surement basis from the desired one as an experimental 
error and deal with it accordingly. 

In the above discussion, we have ignored any decoher- 
ence other than that associated with the measurement- 
induced projection. If the measurement time is much 
smaller than the timescale of decoherence caused by other 
mechanisms, the measurement is completed with mini- 
mal effects of any additional decoherence. Our results 
are valid in this case. A number of different types of de- 
tectors, including QPCs, are approaching this limit where 
the decoherence is limited by the measurement [ll| . in- 
dicating that our results could be observable in such sys- 
tems. 

Quantum state tomography. — One example of a proce- 
dure where the uncontroUability of the measurement ba- 
sis can be harmless is quantum state tomography (QST). 
In QST, one is given a large number of copies of the same 
quantum state, and the goal is to deduce this state. Typ- 
ically, one measures the operator for one third of the 
copies, and similarly for ay and az- Once the average 
values ((7a;), ((7y) and (cTz) are known, the density matrix 
is reconstructed straightforwardly: 

P = ^ (1 + (ctx) 0-x + i^v) o-y + (az) az) ■ (13) 
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In the present case, the measurement bases are cho- 
sen stochastically, and in principle no two of them are 
the same. One must therefore reconstruct the unknown 
quantum state using a procedure that allows for data 
taken from measurements made in any combination of 
bases. One such procedure is the minimization of the 
function 

r(r,0,0) =^[1- r cos (14) 

j 

where r, 6 and are the spherical coordinates of a point 
in the Bloch sphere; j is an index labelling the differ- 
ent runs of the experiment; the direction {Oj,(j)j) defines 
the qubit state obtained in a given measurement; and 
V,{9,(l),9j,(pj) is the angle between the directions (0,(1)) 
and {9j Assuming that the measurement bases cover 
all possible directions (see e.g. Fig. 2c), the convergence 
of this procedure to the correct density matrix should be 
similar to the convergence of the standard QST proce- 
dure explained in the previous paragraph. 

We have simulated QST by repeating the measure- 
ment procedure a large number of times, obtaining a set 
of measurement results (in the form of pre-measurement 
qubit states), and then minimizing the function T(r, 9, </>) 
with the respect to r, 6 and </>. We have chosen several 
initial states covering the Bloch sphere, and the tomogra- 
phy procedure consistently produced the initial state of 
the qubit for the parameters of Figs. 2(b,c). For strong 
qubit-detector coupling [see Fig. 2(a)], the procedure be- 
comes unreliable, because the vast majority of the mea- 
surements are performed in one basis. 



V. CONCLUSION 



In conclusion, we have analyzed the question of what 
information can be extracted from the output signal of 
a QPC that weakly probes the charge state of a charge 
qubit when the qubit Hamiltonian induces oscillations 
between the different charge states. We have shown that 
the measurement basis is determined stochastically every 
time the measurement is repeated. In the case of weak 
qubit-detector coupling, the possible measurement bases 
cover all the possible directions. The measurement basis 
and result can both be extracted from the QPC output 
signal. We have also shown that the information acqui- 
sition rate is independent of the angle f3 between the 
direction defining the charge basis and that defining the 
qubit Hamiltonian. In other words, given enough time, 
the detector will produce a high-fidelity measurement re- 
sult, regardless of the value of (3. These results show 
that, under suitable conditions and by proper analysis, 
the detector's ability to obtain high-fidelity information 
about the state of the qubit is not affected by the compe- 
tition between the measurement and coherent-precession 
dynamics. More detailed analysis of the results discussed 
in this paper is presented elsewhere [T^ . 
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